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1. Introduction 



There are just a few processes, yet to be observed at the LHC, which require a true NNLO 
analysis matching a foreseen experimental precision at the percent level. One of the most notable 
such processes is top-quark pair production. Since the current NLO prediction has a scale depen- 
dence of the order of 13% [jl|] and soft-gluon resummations [j2|, |], ||] are not reliable, the case for 
conducting a complete analysis at the next order of perturbation theory is well grounded. 

Here, I review results, which have been obtained thanks to recent advances in the understand- 
ing of the application of the factorization theorem to collinear limits []|] on the one hand, and to 
progress in the evaluation of massive amplitudes started in [|], f7|] on the other. By combining the 
methods, the leading asymptotics in the high-energy regime of the virtual amplitudes have been 
derived in g]. From then on, the direct evaluation approach has also provided power corrections 
in the mass and fully numerical solutions that cover the whole kinematically interesting range [|l0|]. 



2. Leading asymptotics 

Based on experience, it is relatively clear that the logarithmic terms of the expansion of any 
observable in any limit should be much easier to obtain than the full result. How to actually derive 
them in the case of collinearly singular limits has been first worked out for Bhabha scattering 



[11]. Later on, it turned out that even the constant should be within reach thanks to the use of 



factorization properties of amplitudes [12]. For QCD, the full procedure has been presented in 
[]|] (for application in the context of Bhabha scattering see []l3|]). Until now, the only part that 
poses problems is the one due to contributions of closed heavy particle loops, since they actually 
generate soft singularities, which are process dependent in the sense, that they involve a non-trivial 
color flow between different external legs. 

Let us expand the top-quark pair production amplitude in the strong coupling constant as 
follows 



= 4kcc s 



+ (S)^ (,))+ (i) 2| ^ <2)>+ ^ 



(2.1) 



If the production process involves quark pair annihilation, one can perform the following color 
decomposition of color averaged amplitudes 



2Re(^ (0) |^ (2) ) = (2.2) 
2(N 2 - 1 ) [N 2 A + B + ^ C + Nnfli + Nn h D h + ^E t + ^E h + {n, + n h ) 2 F^j , 

where N is the number of colors, n h the number of heavy quarks and n t the number of light 
quarks. As noted above, all terms which are free of heavy-quark loops, i.e. A,B,C,Di,Ei can be 
obtained from factorization. Although the term F involves heavy quark loops, it is simple enough 
to be inferred analogously to the remaining ones. In principle, to complete the expression, one 
would only have to compute D/, and E%. These contributions are particularly easy to evaluate with 
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the help of methods based on Mellin-Barnes representations [ 14, 15, 16]. The results for all the 
coefficients have been presented in [@], where the authors also checked the factorization formalism 
by computing directly the A, D/ and F coefficients finding, as expected, perfect agreement. The 
interested reader is encouraged to consult [Jsj] for a detailed description of the techniques used. 
Here, I should only stress that the evaluation of the subleading color coefficients is substantially 
more complicated due to the occurrence of non-planar graphs. In the mean time, however, all the 
coefficients have been obtained by the present author, confirming the published results. 

Similarly to the quark annihilation channel, the averaged amplitude in the gluon fusion channel 
can be decomposed as follows 



2Re(^ (0) [.# (2) ) = (N 2 -l)^N 3 A + NB + ^C + ^D + N 2 n,Ei+N 2 n h E h (2.3) 

2 2 \ 

+n,F, + n h F h + ^G, + ^G h + Nn?H x + Nn,n h H lh + Nn h 2 H h + j^I t + -j±I lh + ^4 J ■ 

The results for all the coefficients have been given in The direct computation covered 
AjEhEfaFfaGfaHijHfaHfaliilittJh, whereas the factorization formalism provided A, B,C, D, E\, 
Fi, Gi, Hi, Hih, Hh, I[, Iih, 4. Note that this time it was impossible not to compute non-planar 
graphs, since part of the heavy quark loop contributions involved them. The direct computation 
of the coefficients covered by factorization only is still under way, since it is necessary for the 
evaluation of the full amplitude. 



3. Power corrections 



Due to partem fluxes, a substantial part of top quark pairs will be produced relatively close to 
threshold. It should, therefore, be clear that the high-energy asymptotics described in the previous 
section is insufficient for practical purposes. An exact evaluation of the amplitude seems to be out 
of reach at present, but one can expand in the mass. Introducing 

s={p\+Pi) 2 , t = (pi- p 3 ) 2 -m 2 , x=—, (3.1) 

s 

where p\ and p2 are the momenta of the incoming partons and p 3 is the momentum of the outgoing 
top-quark, and 

xe [1/2(1 -fi), 1/2(1 +j3)], p = ^\-4m 2 /s. (3.2) 

With r = 4mf/s and x = 1/2, which corresponds to 90 degree scattering, the expansion of the 
bare leading color term A is 

36. 4466 + 36. 6376r-2. 11948 r 2 + 0. 318695 r 3 +1.8244 r 4 + 
3.25332 r 5 + 4.54712 r 6 + 5. 73099 r 7 + 6. 82737 r 8 + 7. 85347 r 9 + 8. 82223 r 10 + 
(-7.13933 - 13.2047 r - 1.679 r 2 - 1 . 79336 r 3 - . 682462 r 4 + . 344393 r 5 + 

1.31353 r 6 + 2.23685 r 7 +3.12317r 8 + 3. 97884 r 9 + 4.80856 r 10 ) Log[ms] + 
(4 . 44113 + 2 . 67214 r - 1 . 97207 r 2 - . 387011 r 3 - 0.0835894 r 4 + 0. 0520038 r 5 + 

0.164006 r 6 + . 287566 r 7 + . 428835 r 8 + 0.5859 r 9 + . 755209 r 10 ) Log[ms] 2 + 
(-0 . 608904 - 1 . 10606 r + . 155183 r 2 + 0. 298532 r 3 + 0.425373 r 4 + 0.543219 r 5 + 

0.651184 r 6 + 0.75001 r 7 + 0. 840853 r 8 + 0. 924882 r 9 + 1.00314 r 10 ) Log[ms] 3 + 
(0 . 0625 + . 0833333 r + . 0208333 r 2 + . 00260417 r 3 - . 0078125 r" - 0.0157064 r 5 - 

0.0222168 r 6 - 0.027832 r 7 - 0.032814 r 3 - 0.0373214 r 9 - . 4 1 4 5 8 1 r 10 ) Logfms]" 

(3.3) 
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Figure 1: Bare leading color two-loop amplitude for top quark pair production in the quark annihilation 
channel expanded in the mass. The more divergent terms at threshold (left of the plots) correspond to higher 
orders of expansion. The variables are defined in the text. 



This result in itself looks promising, because we observe asymptotic convergence of the series 
even at threshold. Unfortunately, the subleading color coefficients cannot behave so nicely, as they 
contain Coulomb singularities. Moreover, the series suffers from convergence problems for small 
and large angles as depicted in Fig. |l[ 

Further discussion and details can be found in JToj] . 

4. Numerical solutions 

As the discussion of the previous section shows, it is impossible to proceed further without 
evaluating the amplitude in the threshold region. It turns out that this can be done efficiently with 
the method of differential equations Jl^, |(]] (for previous applications to physical processes 
see [ ^T| , |22|]). The results for the bare leading color coefficient are shown in Fig. ||[ For the details 
the reader is deferred again to [|To|]. 

5. Conclusions 

Although there is still some work to be done before obtaining a true NNLO prediction for 
the total cross section for the top-quark pair production cross section, we are at the stage where 
the virtual amplitudes are almost known (notice that the one-loop squared contributions are known 
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since [23, 24]). The real-virtual corrections are mostly known thanks to the tt+jet calculation [25] 



and the necessary phase space integrations of the real radiation corrections can be done with the 



help of sector decomposition [26, 27], which would then complete the program. Let me also note 



that similar progress has been obtained in the case of gauge-boson pair production [Eq] (here, the 



gauge boson pair production with an additional jet has been studied very recently in [29, BQ]). 



References 

[1] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303 (1988) 607. 

[2] E. Laenen, J. Smith and W. L. van Neerven, Nucl. Phys. B 369 (1992) 543. 

[3] R. Bonciani, S. Catani, M. L. Mangano and P. Nason, Nucl. Phys. B 529 (1998) 424. 

[4] N. Kidonakis and R. Vogt, Phys. Rev. D 68 (2003) 1 14014. 

[5] A. Mitov and S. Moch, JHEP 05 (2007) 001. 

[6] M. Czakon, J. Gluza and T. Riemann, Phys. Rev. D71 (2005) 073009. 
[7] M. Czakon, J. Gluza and T. Riemann, Nucl. Phys. B751 (2006) 1. 
[8] M. Czakon, A. Mitov and S. Moch, Phys. Lett. B651 (2007) 147. 
[9] M. Czakon, A. Mitov, and S. Moch, arXiv:0707.4139 [hep-ph]. 
[10] M. Czakon, arXiv:0803.1400 [hep-ph]. 

[11] E. W. N. Glover, J. B. Tausk and J. J. Van der Bij, Phys. Lett. B 516 (2001) 33. 
[12] A. A. Penin, Phys. Rev. Lett. 95 (2005) 010408. 
[13] T. Becher and K. Melnikov, JHEP 0706 (2007) 084. 



5 



Mass effects and four-particle amplitudes at the two-loop level in QCD 



Michal Czakon 



[14] V.A. Smirnov, Phys. Lett. B460 (1999) 397. 

[15] J.B. Tausk, Phys. Lett. B469 (1999) 225. 

[16] M. Czakon, Comput. Phys. Commun. 175 (2006) 559. 

[17] C. Anastasiou and A. Daleo, JHEP 0610 (2006) 031. 

[18] A. V. Kotikov, Phys. Lett. B 254 (1991) 158. 

[19] E. Remiddi, Nuovo Cim. A 110 (1997) 1435. 

[20] M. Caffo, H. Czyz and E. Remiddi, Nucl. Phys. B 634 (2002) 309. 

[21] R. Boughezal, M. Czakon and T. Schutzmeier, JHEP 0709 (2007) 072. 

[22] M. Czakon and T. Schutzmeier, arXiv:0712.2762 [hep-ph]. 

[23] J.G. Korner, Z. Merebashvili and M. Rogal, Phys. Rev. D73 (2006) 034030. 

[24] J.G. Korner, Z. Merebashvili and M. Rogal, arXiv:0802.0106 [hep-ph]. 

[25] S. Dittmaier, P. Uwer and S. Weinzierl, Phys. Rev. Lett. 98 (2007) 262002. 

[26] T. Binoth and G. Heinrich, Nucl. Phys. B 585 (2000) 741. 

[27] T. Binoth and G. Heinrich, Nucl. Phys. B 693 (2004) 134. 

[28] G. Chachamis, M. Czakon and D. Eiras, arXiv:0802.4028 [hep-ph]. 

[29] S. Dittmaier, S. Kallweit and P. Uwer, Phys. Rev. Lett. 100 (2008) 062003. 

[30] J. M. Campbell, R. Keith Ellis and G. Zanderighi, JHEP 0712 (2007) 056. 



6 



